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Abstract 

We study the large- limit of adjoint fermion one-matrix models. We find one-cut 
solutions of the loop equations for the correlators of these models and show that they 
exhibit third order phase transitions associated with m-th order multi-critical points 
with string susceptibility exponents 7str = — l/w.. We also find critical points which can 
be interpreted as points of first order phase transitions, and we discuss the implications 
of this critical behaviour for the topological expansion of these matrix models. 
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Hermitian matrix models are the classic example of a = quantum field theory where 
'tHooft's topological large- expansion [|I| can be solved explicitly 0, ^. They have recently 
been of interest in the study of the statistical mechanics of random surfaces particularly 
for non-perturbative approaches to lower dimensional string theory 0. There, the large- 
expansion coincides with the genus expansion and the large- iV limit exhibits phase transitions 
which correspond to the continuum limits of the discretized random surface theories 0. 
Unitary matrix models also play a role in 2-dimensional QCD P], mean-field computations 
in lattice gauge theory 0] and various other approaches to higher dimensional gauge theories 
such as induced QCD In this Letter we study a matrix model where the degrees of freedom 
are matrices whose elements are anticommuting Grassmann numbers. The partition function 
is 

Z = Jd^dtP e^' ^^^^^ (1) 
where ^ is a polynomial potential, 

k=i 

i/j and are independent N x N matrices with anticommuting nilpotent elements and the 
integration measure, dip dip = Ylijdipij dipij, is defined using the usual rules for integrating 
Grassmann variables, / dipij ipij = 1, / dipij 1 = 0. We normalize all traces here and in the 
following as V = j^J^iVu- 

Matrix models of this kind have been studied recently by Makeenko and Zarembo , and 
Ambj0rn, Kristjansen and Makeenko |T^. They are motivated by models of induced gauge 



theories using adjoint matter where the Yang-Mills interactions of gluons are induced by loops 
with heavy adjoint scalar fields ||^ or other kinds of matter such as heavy adjoint fermions 
[P, |lT|. Makeenko and Zarembo have shown that the adjoint fermion matrix model (|l]) 
has many of the features of the more familiar Hermitian one-matrix model in the large-A^ 
limit |12] , including multi-critical behaviour with a third order phase transition and string 
susceptibility with critical exponent 7str = — m G !^"^. They also showed that the loop 
equations for the model (P are identical to those for the Hermitian one-matrix model with 
generalized Penner potential [1^ 



l[d(j)ij e-^'*'^(^('^)-2i°g*) (3) 



However, the loop equations for the two models should be solved with different boundary 
conditions and the solution beyond the leading order in large- is different in the two cases. 



It was argued in and |I0|] that, as the Penner model corresponds to a certain statistical 
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theory of triangulated random surfaces, the fermionic matrix model corresponds to a similar 
theory where the genus expansion has alternating signs. The resulting convergence of the 
sum over genera is reflected in the feature of the fermionic matrix model (|I|) that its partition 
function and observables are well-defined since the integrals over Grassmann variables always 
converge, in contrast to Hermitian matrix models |12|] where, typically, in the region of 
interest the integration over Hermitian matrices diverges (reflecting of course the divergence 
of the genus sum in the random surface model). 

In fact, for any polynomial potential the integration over anticommuting variables in (|l]) 
can be formally performed by inserting the matrix- valued delta function 1 = / Hij d(f)ij 6{(f) — 
ipip), where is a Hermitian matrix, and using the identity 



/ TT ^ [d^dij e''^ '^f-^^^ ^ [tt^ det^(-a) e^ = det^ 6{<j>) (4) 

J Ij- 271 J J Ij- 2n \ O0 J 



to obtain 



Z = det^ ( S-] e^' 



0=0 



(5) 



Similarly, any correlator is given by 

det^ (^) tr tr (f)P^ ■ ■ ■ e^' 



< tr i^jj^pY' tr {^^pY^ ■ ■ ■> 



^N^ tr 1/(0) 



(6) 

0=0 



In spite of this good convergence of the partition function, it has been shown that in the infinite 
N limit, the model has third order phase transitions 0]. In this Letter, we shall examine 
the third order transitions in more detail and also argue that there may be a first order phase 
transition. We show that in a simple model with potential VIi/jiIj) = ipip + ^{i/jiIj)^ there are 
two critical points, one of them at the third order phase transition a.t g = Qc = 2/27 which 
was studied in ||, |10| and the other at zero coupling g = 0. In the region g G [0,(7^5 multi- 
branch one-cut solutions exist as well as a multi-cut solution of the matrix model whereas 
when g ^ [0, gd a unique one-cut solution with real free energy exists. 

At g = gc the unique one-cut solution for g > gc connects continuously with the one of the 
3 one-cut solutions which has the minimal free energy in < g < gc- On the other hand, this 
solution is not continuous with the unique one-cut solution for g < 0. The one-cut solution 
for g < connects to one which does not have minimal free energy. Ordinarily, the infinite 
energy barrier (the height of the barrier is of order A^^) between the metastable and stable 
one-cut solutions prevents tunneling and also a phase transition from occuring. However, 
if we restrict attention to one-cut solutions, and follow them over the range of g, we must 
encounter a discontinuity of the free energy somewhere, i.e. a first order phase transition. 
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The choice of stable solution (the one-cut solution with minimum free energy) in the regime 
g G [0, results in a first order phase transition aX g = {]. With this stability requirement 
perturbation theory near the Gaussian point g = Q does not correctly reflect the properties of 
the theory when g is small and negative. The critical point gc corresponds to the usual m = 2 
multi-critical point with 7str = ~l/2. In the following, we shall also explicitly construct one- 
cut solutions for generic symmetric potentials and thus show how higher order multi-critical 
points can be realized in these models. 

First, we shall discuss some general properties of the adjoint fermion matrix model (HI). 
It possesses a continuous symmetry ip UipV~^, ip VipU~^ with {U,V} G GL{N,€) ® 
GL{N,(C). In spite of this large degree of symmetry, it is not possible to diagonalize a matrix 
with anticommuting entries. Thus, unlike the more familiar Hermitian one-matrix models, the 
model (|1|) cannot be written as a statistical theory of eigenvalues. Nevertheless, in the large- iV 
limit it shares many of the properties of such a theory |^ ^. Furthermore the large degree of 
symmetry restricts the observables to those which are essentially invariant functions of ipip. 
The chiral transformation ip ^ ip, ip ^ tr [ipip)^ (—1)^+^ tr [ipip)^ is a symmetry 

when the potential is an odd polynomial. Furthermore, in that case, all even moments vanish, 

< tr [ipilj)^'' > =0. This is the analog of a symmetric potential for a Hermitian matrix 
model If. 

In the large- limit, correlators of the matrix model factorize 

< tr fi^tp) tr g{^^) > = < tr fiipip) > < tr gitp^) > +0{1/N'^) (7) 

This factorization property follows from the existence of a finite large- limit for the corre- 
lators < tr [ipip)^ > for arbitrary polynomial potential VIi/jiIj) = J2k ^{'4'i^)^i since then the 
connected correlators are given by 

< tr (#)P tr {,p,ljf —p— < tr {^P^^f > ^ j;p (§) 

Factorization and symmetry imply that the large- limit of the model is completely char- 
acterized by the set of correlators < tr [ipip)^ >. Since the signs of the actions in (P and 
(1^) are opposite (from the different boundary conditions), it also follows that the connected 
correlators of the fermionic matrix model alternate in sign relative to those of the generalized 
Hermitian Penner model (^. It is this property that leads to an alternating series for the 
large- iV genus expansion of the fermionic matrix model ([l|). 

When is finite, nilpotency of the components of ip and ip implies that the moments 

< tr [ipipY > are non-zero only for k < N"^ and are therefore finite in number. In this case. 
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the moment generating function 



1 ^ 1 



u{z) = < tr > = E < tr i^^)' > ^ (9) 

has singularities only at the origin in the complex 2;-plane. The finite set of moments can 
always be obtained from a (not unique) distribution function p with support in the complex 
plane 

< tr (V^V)^ > = J da p{a)a'' with J da p{a) = 1 (10) 

The support of p can be deduced from the position of the singularities of in (P). When N 
and therefore the number of moments is finite the support of p is concentrated near the origin 

In the large- limit, the spectral function p{a) can be a function with support on some 
contour in the complex plane. The distribution function p is the analog in the fermionic 
matrix model of the density of eigenvalues in Hermitian one-matrix models as the quantity 
which specifies the solution of the model in the infinite N limit 

The generating functions for the connected correlators are 

11 

UJn[Zl, ...,Zn)= < tr ^ • • ■ tr ^ >conn (12) 

Zl -Iplp Zn- WW 

In particular, the loop correlator 

ujAz) =oj{z) = I da (13) 
J z — a 

is analytic in z away from the support of p in the complex plane. The distribution function 
can be determined by computing the discontinuity of uj{z) across its support, where 

uj{\ ± e^) = -Ida f =F T-^7rp(A) for A G supp p (14) 
J A — a \e±_\ 

where e^(A) is a complex number with infinitesimal amplitude and direction perpendicular to 
the integration contour at the point A. The Schwinger- Dyson equation for the loop correlator 
uj{z) can be derived from the invariance of the partition function (0) under arbitrary changes 
of variables. It can be cast in the form P, T(]| 



- d 
dip dijj 



i/jijj 



(15) 
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In contrast to Hermitian matrix models |T2[ , the identity ( ]T5| ) is exact for fermionic matrices. 
It leads to 

- zuj{zf + (2 - zV'{z)) uj{z) + V'{z) + P{z) = zuj2{z, z) (16) 
where Pi^z^ is a polynomial of degree K — 2 

P{z) = j^gu J2 < i^^)"-'-" > (17) 

k=2 p=0 

Factorization implies that the connected correlators are all suppressed by factors of 1/A^^ 
and the term on the right-hand side of the loop equation (|16D vanishes in the large- limit. 
Then the loop equation has solution 



/ , 1 y\z) 1 



' zV'{z)\^ 



+ zP{z) 



where the sign of the square root is chosen to yield the correct asymptotic behavior lo{^z) ^ 1/ z 
at \z\ oo. The branches of the square root must be placed so that it is negative near the 
origin in order to cancel the pole at z = 0. If the potential is a polynomial of order the 
solution will in general possess a square root singularity with K branch cuts and the spectral 
density p will have K contours in its support. 

The simplest solution of the model is the one-cut solution which assumes that the singu- 
larities of uj{z) consist of only a single square root branch cut, so that the distribution p has 
support only on one arc in the complex plane with endpoints at some complex values ai and 
a2. This solution for uj{z) can be represented in the form |P, |10| 



. dw V'iw) - 21 w 

uj{z) 



Tc Am z — w \ {w — ai){w — 

where the closed contour C encloses the support of the spectral function but not the point 
w = z. The endpoints of the cut can then be found by imposing the asymptotic boundary 
condition uj{z) ^ 1/ z ai 1^1 — oo on the solution (|19]), which leads to the two equations 

dw V'{w) -2/w ^ / dw wV'{w) -2 



^ V (""^ ~ ai)(w — 02) ^ yi'^~ ai){w — 02) 

To determine the precise location of the support contour of p in the complex plane, we first 
use the observation of @ that the large-iV equation ( |16|) for the loop correlator is identical 
to the loop equation for the generalized Penner model (|^). In the large- limit, the spectral 
density therefore obeys the saddle-point equation 

2/X-V'{X) r p{a) 



-Ida , A G supp p (21) 

J X — a 
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Note that this equation can be obtained from the discontinuity (|I4|) of the loop correlator 
([T8|), and it also follows from the minimization condition for the free energy 

F = lim log Zp= f da pia) (Via) - 2 log a) + Uda dX p(a)p(X) \og(a - A) (22) 

Af^oo A* ^ J J J 

with respect to the distribution function p. Note the change in sign of the fermionic free 
energy relative to the Hermitian case. The double integral in (|2^) is evaluated by integrating 
up the saddle-point equation (pl|). This introduces a logarithmic divergence at A = arising 
from the Penner potential in (^ which we remove by subtracting from (|2^ ) the Gaussian free 
energy Fq defined by setting (yf^ = for > 1 in (^2]) 



F-Fg = ^J da p{a) {V{a) - 2 log a) + -J^da p{a) log a (23) 

where we have ignored terms independent oi Qk, k > 1. The support contour of p can be 
determined from the David primitive function |n 

G{w) = £ dz - V'{z) - 2uo{z)^ (24) 

The support of p is an arc connecting ai to 02 in the complex plane along which G{w) is 
purely imaginary and which can be embedded in a region where Re G{w) < [0 . 

For illustration, we shall consider the cubic potential 



V{z) =tz + ^z^ (25) 



for which 



I t gz"^ 1 



uiz) = + — V^?'^' + 2t^?^^ + {t^ + M)^' + 4 (26) 

where ^ is the as yet unknown correlator ^ = < tr ipil) >. In this case the vanishing of all 
even moments, / da p{a)a'^^ = 0, implies that the endpoints of the support contour of the 
continuous function p lie in the complex plane and are symmetric on reflection through the 
origin. Furthermore, an application of Wick's theorem shows that the series in the odd 
moments is alternating. 

Generically the square root in uj{z) has three branch cuts, so that in the general case the 
distribution p will have three disjoint and symmetric (about the origin) support contours. 
The one-cut solution for (p6D takes the form 

^(^) = ~~2~^ 2z V ^ ^ ^ 



where comparing the polynomial coefficients in (pTl) with those of (^) shows that the param- 
eter h and the correlator are determined by the two equations 

±h^ -tl? + 2g = Q (28) 
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- {t"^ + 4:g^)b ± 8g = (29) 

The sign ambiguity here can be ehminated by requiring that ai g = the correct Gaussian 
value b{g = 0,t) = t for b be attainable. This is the boundary condition that is relevant 
for an interpretation of this matrix model as a discretized random surface theory, i.e. for a 
consistent perturbative expansion of the model in the coupling constant g. It means that we 
take the positive sign in the above equations. The choice of negative sign yields solutions 
with boundary conditions aX g = appropriate to generalized Penner models [0, |13[. The 
equation (p8|) and this sign ambiguity also follow from the contour integrals (|20|). 

We assume henceforth that t is a positive constant. The 3 solutions of (|28|) are 

b,{x,t)=^-{(3'/%x)+(3-'/%x) + l) (30) 
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where 

P{x) = 2x-l + 2^x{x - 1) (32) 

and we have introduced the scaling parameter x = 1 — . When x < {g > gc = ^) 
or X > 1 {g < 0), (3{x) is a. monotone real- valued function with (3{x) > 1 for x > 1 and 
/3(x) < — 1 for X < 0. In the region < x < 1 {0 < g < gc), P{x) is a complex- valued function 
with unit modulus. The function (^) is always real-valued and the region < x < 1 is the 
region wherein all 3 roots (^), (|3l|) of the cubic equation ( p8l) are real. These 3 roots can all 
be obtained from (|30D by choosing the 3 inequivalent cube roots of /5(x). For x ^ (0, 1) the 
solutions (^) are complex. For the fermionic matrix model, where the distribution function p 
can be complex-valued, there is no immediate reason to disregard complex-valued endpoints 
for the support of p. However, the free energy ( ^3]) for the cubic potential (|25|) up to terms 
independent of b and g is 

where we have used the spectral density determined by (jM]) and ( p7|) 



P(a) = ^(6 + ^7«^)^1 + ^ ; aea (34) 

where Cb = supp p and b{g, t) are given by (0) and (|3TD . It is immediately seen that the free 
energy (p3D is complex- valued for the values (|3ll) of b{x,t) for x ^ (0, 1). Such a free energy 



leads to an unstable state and we therefore consider only the real- valued solutions to (p8|). 
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The support contour Cf, on which 
which is 



G{z) = -\/W-z^ + 4 - sgn (6) log 



is defined is found from the David function ( p^ 




3|6|3 



6V + 4f^' 



zvr sgn h 



(35) 



where the branch of the square roots in ( ^5|) is taken to be the straight fine joining the 
points ±2i/|6|. A careful examination of the equation Re Gi^z) = and of the region where 
Re G{z) < shows that the contour Cb cannot cross the imaginary axis for | Im 2;| > 2/\b\ and 
that it crosses the real axis at some non-zero values of order ±1/|6|. The regions Re G{z) < 
are to the right of these crossing points (but note that Re G{z) changes sign across Cb). Thus 
the contour Cb in (|3^) can be taken as the counterclockwise oriented half-circle of radius 
2/\b\ in the first and fourth quadrants of the complex a-plane. It is easy to verify that with 
this definition of p the equations (0) and (|T^) are satisfied, as is ( pQ]) from evaluating the 
correlator ^ = J da p{a)a with this distribution function. 

There are 2 critical points in this large-A^ matrix model, ai g = and g = gc, which 



separate 3 phases determined by the analytic structure of the function (|32D , i.e. the one-cut 
solution is a non-analytic function of x about x = and x = 1 where it acquires a square root 
branch cut. For x > 1 the solution 



bo{x,t) 



1/3 



1+ (2x-l + 2^Jx{x-l)] + (2x - 1 + 2^a;(a; - 1) 



-1/3 



X > 1 (36) 



of ( P^D satisfies the Gaussian boundary condition bo{x = l,t) = t. When x < the real 
solution for b is 



bo{x,t) 



1 - 



2x-l + 2Jx{x-l) 



1/3 



2x - 1 + Jx(x - 1 



-1/3- 



X < 



(37) 



As X is varied between and 1, (3{x) has modulus one and phase which varies from vr to 
0: P{x) = e^'^^^^ where 

/ 9. /-T-n — t\\ 

(38) 



2Jx(l -x) 

(x) = arctan | — G [0, vr] 

2x — 1 



The arctangent function in (|38|) is well-defined only up to an integral multiple of 27r, and the 
three real solutions for b are 



bix,t) = - 



1 /2Jx(l-x) 

1 + 2 cos < - arctan — 

3 I 2x - 1 



2mT 



0<x<l, n = 0,1,2 (39) 
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The branch which matches ( |37D is the one with n = 1, whereas the branch which matches 
( p6| ) is the one with n = 0. The branch with n = 2 does not connect with either solution. 
Any of these 3 branches can be used to define the one-cut solution (pT]). The free energy ( |33D 
is positive for all x G (0, 1) for the n = branch, negative for all x G (0, 1) for the n = 1 
branch, and for the n = 2 branch it is positive for < x < | and fiips sign for the rest of 
the interval at x = |. The n = 1 branch in ( p9| ) is the ground state solution in the region 
< X < 1. 

The free energy associated with this stable one-cut solution is discontinuous across g = 0, 
and thus with this choice of branch in the regime < x < 1 the Gaussian point of this matrix 
model is a critical point of a first order phase transition. The other one-cut solution which is 
a perturbation of the Gaussian solution is metastable but can still be thought of as a valid 
solution of the model since the barrier is infinite at = oo. This is similar to the situation in 



the Hermitian one-matrix model with symmetric polynomial potential of degree 6 . There 
a phase transition occurs due to an infinite volume effect, as opposed to a large- effect 
where the only possibilities could be second or third order phase transitions. There is also the 
possibility that the loop correlator (^) evolves into a three-cut phase at g = 0, corresponding 
to a third order phase transition at g = (see below), but there is no immediate indication 
of this since in the fermionic case the spectral measure p{a) da need not be positive. This 
possibility is also suggested by the exact form (^61) of the loop correlator: Although the one- 
cut ansatz ( pTl) is insensitive to a change in sign of g, the analytic properties of (|2^) might be 
affected by the passage through g = 0. 

The existence of 3 phases in this matrix model and the possibility of a first order phase 
transition at the Gaussian point g = are completely unlike what occurs in the conventional 



polynomial Hermitian matrix models [§, O or in Penner models [10, fq|. The unusual phase 



transition aX g = would imply that the consistent perturbative expansion of the theory in g 
has a preferred direction through values of g with a definite sign (corresponding to — sgn t). 
This fact is important for the interpretation of the fermionic one-matrix model as a statistical 
theory of discretized random surfaces. Notice, however, that the free energy (^) with the 
choice of stable branch for x G (0, 1) is continuous across the point g = gc which coincides 
with the critical point found in ||10|| . 

The scaling behaviour of the matrix model in the vicinity of its critical points is determined 
by the string susceptibihty 

The critical exponent 7str each critical point (7^*^ is defined by the leading non-analytic 
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behaviour of (|0|) § 



X{9,t)r^s{g-gi'^r''^^^ as g^gi^^ (41) 



where denotes the most singular part of the function in a neighbourhood of the critical 
point. In terms of the scaling variable x, the susceptibility (j^) is 



/ , 972 



t6(/52/3 + /5I/3 + 1)6 



1 - 8a; + 8a;^ + 4./—^ (l - 3s + 2^^ 
V a; — 1 



(42) 



From (^) we find that the leading singular parts of the susceptibility near each of the two 
critical points g = gc and g = are respectively 

X{x,t) r^s - ^^^q'^ Vx as x^O (43) 
11648 



X{x,t) ~5 -^^Vx - 1 as x-^1 (44) 
Both critical points therefore have string constant 7str = ~l/2 which coincides with those of 



the usual genus zero m = 2 quantum gravity models ||T2 . 

In particular, this shows that with the choice of stable branch for x G (0, 1) the phase 
transition at the non-zero critical coupling g = gc is of third order. Indeed, the critical point 



g = gc enjoys all of the properties of a conventional m = 2 multi-critical point |T^. The 



information relevant to the topological genus expansion in string theory is contained in this 



non-zero third order critical coupling. As argued in |]I0|, this m = 2 multi-critical point 
results in a genus expansion which is an alternating series but otherwise coincides with the 
usual Painleve expansion 0. 

Thus the 2 critical points of the fermionic one-matrix model, which arise as those points in 
parameter space where the cubic equation (^) which determines the one-cut solution acquires 
a double real root, provide the necessary information for the perturbative and topological 
expansions of the theory. In the 2 phases outside of the region < a; < 1 the one-cut solution 
is unique, while in the phase x G (0, 1) a three-cut solution can in addition exist. The scaling 
behaviours (|43|) and (|^) indicate that the 2 transitions into the multi-cut phase would both 
be of third order, while the transitions into the stable one-cut phase are of first and third 
order. The existence of a single-cut or multi-cut phase in < a; < 1 is determined by which 
one of these 2 possibilities is in fact the vacuum state. The phase structure of this model 
is reminiscent of the triple point of a liquid-vapour-solid system, where there is a classical 
forbidden region in the liquid phase for the isotherms to pass through. The system can 
however supercool from the vapour to solid phase by tunneling through this region which 
corresponds to a first order phase transition, while the other phase transitions are of third 
order. In our case this "forbidden" region is < x < 1 (as this is where the function f3{x) 
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becomes non-analytic) and it has width g^. = 2t^/27. The triple point where all 3 phases 
coexist is then obtained by letting t —>■ 0. 

We conclude by discussing the critical behaviour associated with higher order potentials. 
For simplicity we consider the chirally symmetric case where the potential (0) is a generic odd 
polynomial, i.e. g2k = for all k in (|^), with K = deg \^ > 3 an odd integer. The solution for 
the loop correlator is then 

1/2 



N -L 92k-l 2(k-l) . /I I „ „ ^2(fc+m-l) i „ c 



2mZ 



2k 



V 



k,m=l 



k,m=l 



where ^2m are the as yet unknown moments ^2m 
for p5|) takes the form 

/ 

-L .92fc-l 9l'^_1^ . J- 

UJ(Z 



(45) 

< tr [ipipY"^^-^ >. The one-cut solution 



1 _ 92k-l 2{k-l) , ^ 

- ^ 2 +2^ 



A:=l 



_ftr-i 



+ E «2fc^'' + & + 4/62 (46) 



V 



k=l 



where we have fixed the sign in front of the parameter b by the same convention as before. 
The one-cut ansatz (|i6| ) along with the general solution (^5|) together involve K — 1 unknown 
parameters - b, the {K — 3)/2 polynomial coefficients a2k, and the {K — l)/2 correlators C,2k- 



These parameters can be found by comparing the various polynomial coefficients of (^5|) 
with those of (^61) , which leads to the set of equations 



b^+8a2-b 



2 



gl + E 92k+i^2k \=0 , 2b^+8ba^+Aal-b^ 2g,g, + E ^^2^+36^ = (47) 



V 



k=l 



V 



k=l 



K-2 K-1 

2bgKaK-5 + ^a,K-3-b E 92k-i92K^2k-z = , 2b'^ gKax-^ + ^R-b'^ E 92k-i92K~2k-i = 

k=l k=l 

(48) 



K~l 
2 



K+1 

2 



2b'^aK-3 + 8gK-b gxi + E 92k-i9K- 



2k 



k=l 



, 2b^gK + 8bgKa2-b'^ E 92k-i9K-2k+2 = 

(49) 



k=l 



When K > 7 we have in addition the sets of equations 



m-2 



26^a2(m-l) + 8ba2m + Aa20'2{m-l) + E (^^'^2A:«2(m-fe-l) + ^a2k(^2(m-k) 

\ 



k=l 



K + l-2m 
2 



I E 92k-l92{m-k)+l + E 92{rn+k)-1^2k 
k=l k=l 



(50) 
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for 3 < m < ^^2^, and 

m-2 

2b gKCi2m-K-l + 86a2m-X+l '22A:fl2{m-fe-l) 

k=l 

m—1 m 

+4 ^ a2ka2{m-k) -b'^Yl 92k-ig2(m-k)+l = (51) 

k=l k=l 

for ^ < m < K - 3. 

(^71) -(PH) yield a complete set of equations for the K — 1 unknown coefficients of the one- 
cut solution ( ^61) in terms of the coupling constants of the potential (§). The parameter b can 
alternatively be found from the contour integrals ( pOf ) which lead to a i^-th order equation 
for b 

K+l 

, ± (zl)!g|^,,,_.,.-=.« . (52) 

Since K is odd this equation always has a real solution, and as before the one-cut solution 
can always be constructed. The spectral density is 



K-3 
2 



v 



k=l 




Qkoi^-^ + «2fca^'' 1^/1 + ^ ; « e a (53) 



The first {K — l)/2 moments of this distribution are given by the solutions to (^^-(^Tj). 

In general (|52|) will acquire multiple real roots at some coupling 5'2fc-i which will be a 
critical point of a third order phase transition with string constant 7str = —1/2. We expect 
that the phase with multiple real roots will be bounded by other coupling constant values so 
that the model will contain several critical points corresponding possibly to different order 
phase transitions. These critical points might provide information about this matrix model 
other than that relevant to the genus expansion. Notice that since the potential now depends 
on more parameters, we can adjust them in such a way that m—1 regular zeroes of ( |53D 
coalesce with a cut end-point ±.2i/\b\ at criticality for the same critical coupling g2k~v 92k-i 
will then be an m-th order multi-critical point with susceptibility exponent 7str = —1/m 
[|12| . In this way the construction above shows explicitly how more general two-dimensional 
quantum gravity theories |T^] may be obtained from continuum limits of fermionic matrix 
models. 
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